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Abstract 

The doubly charged scalar boson (H ) is introduced in several models of the new physics beyond the 
Q-[ standard model. The has Yukawa interactions with two left-handed charged leptons or two right- 

i <""] handed charged leptons depending on the models. We study kinematical properties of H ±±: decay products 

through tau leptons in order to discriminate the chiral structures of the new Yukawa interaction. The 
chirality of tau leptons can be measured by the energy distributions of the tau decay products, and thus 
the chiral structure of the new Yukawa interaction can be traced in the invariant-mass distributions of the 
i? ±=t decay products. We perform simulation studies for the typical decay patterns of the ±:t with simple 
event selections and tau-tagging procedures, and show that the chiral structure of the Yukawa interactions 
of 7J ±± can be distinguished by measuring the invariant-mass distributions. 
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I. INTRODUCTION 

The existence of the neutrino masses has been established well [l-6]. However, neutrinos are 
massless in the standard model (SM) because of the absence of the right-handed partners. If the 
lepton number conservation is violated in a new physics model beyond the SM, neutrinos can be 
Major ana particles which form a mass term with its self-conjugation field only [7[, since neutrinos 
are electrically neutral unlike to all other SM fermions. Therefore, it seems natural to expect that 
the possible Majorana nature of neutrinos provides the reason why neutrinos have very different 
masses from those of other SM fermions. 

The doubly charged scalar boson H , which has a twice electrical charge of the electron, exists 
in several models to generate Majorana neutrino masses. For instance, the particle is a member 
of an SU(2)x triplet scalar field in the Higgs triplet model (HTM) Qj. The triplet field develops 
a tiny vacuum expectation value (VEV), which breaks the lepton number conservation and is the 
source of the neutrino mass. Such a triplet field appears also in some models of extended gauge 
symmetries {{J. On the other hand, a doubly charged scalar boson is introduced as an SU(2)l 
singlet scalar field in the Zee-Babu model (ZBM) 10| which generates Majorana neutrino masses 
at the two-loop level. In these models with H , its Yukawa interactions with charged leptons 
depend on the SU(2)x property of H . Namely, H from an SU(2)^ triplet field couples only 
with left-handed charged leptons ij^ while the one from an SU(2)^ singlet field interacts only 



with right-handec 
multiplet scalars 



charged leptons £ R . Furthermore, H can be a component of other SU(2)x, 



111 ] . and such H also has Yukawa interactions with two left-handed or two 



right-handed charged leptons through the mixings between leptons and new fermions. In any case, 
both of two charged leptons which couple with H via the Yukawa interaction are left-handed 
or right-handed. The discrimination of the chiral structure of the Yukawa interaction plays an 
important role to distinguish these models. 

The can be produced by the pair creation process, pp — > 7*/Z* — > H ++ H~~ . For SU(2)x, 
non-singlet representations, the associated production — > W T * — > H^H^ with a singly charged 
scalar boson (H^) is also possible 121. Theoretical studies for decaying into same-signed 

nn 

leptons and weak gauge bosons can be found in, e.g., Refs. [131 . 1141 ]. The experimental search 



results for iJ ±=l= have been available, where purely leptonic decay channels are assumed 15H17]. 
We comment that these bounds on the mass are dependent on the production mechanism, 

the decay branching ratios, and the mass spectrum of the scalar boson multiplets [18 1. 



In this letter, we study the consequence of the chiral structure of the Yukawa interaction (of the 



2 



doubly charged scalar boson with two charged leptons) to the kinematical distribution involving 
the decay of tau leptons. The polarization of r leptons is known to be probed by its decay products, 
and can be exploited to test the structure of new interactions in the models beyond the SM 
In Section II, models of neutrino masses with are introduced with particular attention to the 
chiral structure of the Yukawa interaction. In Section III, the polarization dependences of the 
decay distributions of r leptons are reviewed, and the invariant-mass distributions of final-state 
particles in the decay of ff ±=t into at least one r lepton are discussed. Simulation results including 
r-tagging and simple kinematical cuts are also presented. Conclusions are given in Section IV. 



II. MODELS WITH DOUBLY CHARGED SCALAR BOSONS 

In this section, we briefly present examples of models which include the doubly charged scalar 
boson. 



The first example is the HTM 



fl. 



In this model, an SU(2)^ adjoint scalar field A with hy- 
percharge Y = 1 is introduced in order to generate masses of neutrinos via the triplet Yukawa 
interaction. The new Yukawa interaction is given by 

^htm™ = -L c h M icT2^L + Il.c., (1) 

where L = (u L , t~^) T is the lepton doublet field, the Yukawa coupling matrix is symmetric h M = tij^, 
(Ti(i = 1-3) are the Pauli matrices, and 

/a+A/2 A++ \ 
A= 7 • (2) 

\ A -A+/V2J 

In the HTM, the doubly charged scalar boson interacts with a pair of left-handed charged leptons. 
The neutrino mass matrix in the flavor basis is obtained as M^ 1 ™ = 2/ij w -(A°) = C/mns M u U^ ns , 
where (A ) is the VEV of the triplet field, M v is the neutrino mass matrix in the diagonal basis, and 
Um~ns is the Maki-Nakagawa-Sakata (MNS) matrix for the lepton flavor mixing. Since the neutrino 
mass matrix is directly related to the Yukawa matrix, the decay patterns of the doubly charged 
scalar boson are constrained by observed neutrino oscillation data. For example, {J^m) ^ (^Af ) rr 
and (h M ) e n ~ {h M )er are required because the observed neutrino mass matrix approximately has 
the //-T exchange symmetry. By assuming the realistic values of decay branching ratios, constraints 
on the mass of A ±=l= are obtained > 400GevQ. 

The next example is the ZBM [l(3|. Two SU(2)i singlet scalar bosons, k~ (Y = —1) and 
k (Y = —2), are introduced in the ZBM to generate tiny neutrino masses at the two-loop level. 



The new interaction terms which relevant to the radiative neutrino mass are 



-ZBM 



-L c Y a ia 2 Lk + - (£ R ) C Y s C R k ++ - fx ATAT A; ++ + H.c, (3) 



where Y a = —Yj and Y s = Yj . The doubly charged scalar boson in this model interacts with 
right-handed charged leptons. If a lepton number 2 is assigned to k~ and k , a coupling constant 
/i is the soft breaking parameter of the lepton number conservation. The neutrino mass matrix is 
calculated^ as (M^ bm ) q ^ = 16/i (Y*) a £ m £ (Y s ) U i I U i m £ , (Yj)^g, where the loop function I U i is given 



in Ref. 



22|. In order to describe the observed neutrino oscillation parameters, (^) /i ^(m Al /m T ) 



(Y s ) tlT (m l j_/m T ) ~ (Y S ) TT is favored in the ZBM 



23J. This may suggest that k — > t r t r would 



be highly suppressed while k — > fi R t r could be sizable in the ZBM. Assuming purely muonic 
decay mode, the mass of k ±i: is constrained to be m k ±± > 250 GeV [l^. 1 



III. TAU POLARIZATIONS AND DOUBLY CHARGED SCALAR BOSON DECAYS 
A. Decay distributions of polarized tau leptons 

In this section, we review the polarization dependence of decays of r's, and discuss how that 
could be traced in the case of decay through r's. In the following discussion, we assume 

that the leptonic decays of doubly charged scalar bosons occur via the Yukawa interactions, e.g., 
H^-~£ x (t x ) c (X = L,R), where H~j^~ (H R ~) denotes H only with the left-handed (right-handed) 
interaction. Hereafter, t denotes e or /x. 

First, let us consider the lepton flavor violating (LFV) decay H — > 1~t~ followed by r~ — > 
ix~v. The branching ratio of the pionic decay of r is about 11 % while the branching ratio of the 
total hadronic decay is about 65 %. The invariant-mass of £ir is expressed as Mf n = z m 2 H ±± , where 
m H±± is the mass of and z = E W /E T ; the and E T are energies of a pion and a r lepton 
in the laboratory frame, respectively. 

The distributions of the pion energy fraction z (namely, of the invariant-mass Mf ) are given 

as 

Vl(z) = Fl(z) = 2(1 - z), (4a) 
V R (z) = F%(z) = 2z, (4b) 

If we use theoretical curves in Fig. 2 of Ref. with the result of Ref. [3 (m A ±± 

> 391 GeV for pair-produced 
A with the 100 % decay branching ratio into a muon pair.), we would naively arrive at m k ±± > 320 GeV. 
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FIG. 1: Distributions of the invariant- mass of in (left panel) and of U T (right panel) from the LFV decay 
H — > £~t~ , where £( r ) = e, /i. The invariant-mass distributions through the decay of tl (tr) are plotted 



in the dotted (dashed) curves. 



a 



where F£ R (z) are the fragmentation functions of r L R — > it v decay in the collinear limit |1S||. The 
fragmentation functions for the other hadronic decay modes are also known but less sensitive to 
the polarization of r [l9|. We will utilize these decay modes in the simulation study later. 

When the LFV decay is followed by the leptonic decays of r's, the dilepton invariant-mass is 



express as M] 



z m 



2 H ±±t where z is the energy fraction of the daughter lepton to the parent r 



lepton. We denote £^ as £^ from decays of r^. 2 The total branching ratio of the leptonic decays 
is about 35%. The distributions of z (namely, of Mf t ) are given as 



V{{z) = F[(z) = -(1 - z% 

V R (z) = F R {z) = 2(1 -zf (1 + 2z), 



(5a) 
(5b) 



where R (z) are the fragmentation functions of t£ r — > t~vv decay in the collinear limit 

In Fig. [U we plot the distributions of the invariant-masses of £ir (left panel) and ll T (right panel) 

for the LFV decay H — > £~t~ followed by the pionic and leptonic decays of the r, respectively. 

The invariant-mass distributions via the decay of tl (tr) are plotted in the dotted (dashed) curves. 

For the pionic decay channel, the distributions are linear in z and have an opposite behavior between 

tl and tr. On the other hand, the right panel of Fig. Q] shows that the distribution of for the 

leptonic r decay would be less sensitive to the r polarization than that of the pionic channel. 
Next, we consider the decay mode H — > t~t~ . The decay pattern of the two r's can be 

classified into three categories: hadronic channels (e.g. vr7r), semi- leptonic channels (e.g. £ T ir), and 



The notation tr T \tr T s indicates not only ee and fifj, but also e\x. 
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FIG. 2: Distributions of the invariant- mass of rnr (left panel), £ T ir (middle panel) and of l r l r (right panel) 
from H — > t~t~, where t T ~ e, fj,. The invariant-mass distributions through the decay of tl (tr) axe 
plotted in the dotted (dashed) curves. 

purely leptonic channels (£ T i T ). The distributions of tttt invariant-mass M in H — > t~t~ —> 
■k~it~vi> decay chain are calculated by convoluting the fragmentation functions of the pionic decays 
of r's in Eqs. ((4|) as follows: 



m(*)= [ 1 ^Fl(z 1 )FE(z/z 1 )=4 

J z z l 



(1 + z) log - + 2z - 2 



Vbr{z)= [ 1 ^F^z 1 )F^(z/z 1 ) = 4zlog- 

J z z l z 



(6a) 
(6b) 



where z = M^ n /m 2 H±± and VH (V]? R ) is the distribution for (H^). 

The distributions of the £ t tt invariant-mass M^ n for the H — > t~t~ — 
are given by 

6 log z 6 + 9z 

z 



VU*)= I d -^Fl{z l )Fl{z/z l ) = A - 

Jz z l y 

©feW = f ^ Fh(*i) FR = 4(1 - *) 

J z z l 



? T 7r vvv decay chain 



(7a) 
(7b) 



where z = Ml w /m 2 H±± . The distribution Vf L {V%) is for (H^ 

The dilepton invariant-mass distributions for the H — > t~t~ — > t~t~vvvv decay chain are 
given by 



f dZl n n 16 

where z = M}^Jm\ ±± and V% L {V%) is for {H^ 



2-2z 3 -3(l + z 3 )log-^ 



5 - 27z 2 + 32z 3 + 3(1 + 9z 2 + Az 6 ) log — 

z 



(8a) 
(8b) 
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In Fig. [2 we plot the distributions of the invariant-mass of mr (left panel), £ t tt (middle panel), 
and l r t T (right panel) from H — > t~t~ . The invariant-mass distributions through the decay of 
T L (tr) are plotted in the dotted (dashed) curves. Owing to the large polarization dependence of 
the pionic decay fragmentation function in Eqs. ©, the tttt invariant-mass distribution has a large 
power for the spin analysis. Namely, the distribution takes a large value in the region of a small z 
for the decay while it is large for z ~ 0.4 in the decay. The £ t tt and l T i T distributions 
in the middle and right panels of Fig. [21 respectively, take a large value for a small z for both 
and Ht^, thus are not useful as the r polarization discriminator. 

To summarize, the invariant-mass distributions of Irr (and ££ T ) in the H — > £~t~ decay 
and tttt in the H~ — > t~t~ decay are good for discriminating the r polarization and thus the 
chiral structure of the Yukawa interaction with H^. Notice that the lit (l€) final-state can be 
reached through the other decay chain(s) such like H — > t~t~ — > i~ix~vvv (H — > t~t~ and 
H — > t~t~ — > i~£~vvvv). If these decay modes for exist simultaneously, the signatures 

from these decay chains would mix in general. It should be possible to treat the mixed signatures or 
to divide them by kinematical cuts. However, since such analyses depend on the detail branching 
ratio of the decay, we will not consider them in this study. 



B. Simulation results 



In order to perform realistic studies for collider experiments, we examine a Monte-Carlo simula- 
tion for the H ++ H pair production and their decays up to parton-showering and hadronizations. 
We generate signal events of pp — > H ++ H by using Pythia 



TAUOLA 



241 ] with handling the r decay by 



25( | incorporating the chiral properties of the Yukawa interactions of i? ±± . We consider 



only the leptonic decay of -ff ±=t , and pick up several patterns of the pair of decays suited for the r 
polarization measurement. The m H±± is set to be 400 GeV, and the collider energy to y/s = 14TeV. 
For the analysis, we use lighter leptons (e and /u) with p^ T > 15 GeV and \r\i\ < 2.5, where p T is 
the transverse momentum and 77 is the pseudo rapidity. To find the hadronically decaying r's, we 
perform r-tagging for every jets with plj, > 25 GeV and \r]j\ < 2.5 which are constructed by the 



anti-/ET algorithm 



26| | with R = 0.4. For the r-tagging, we use two methods. The first method 



is devoted to extract the r — > txv decay. Namely, it is tagged as the pionic r-jet (7r r ) if a jet has 
only 1 charged hadron and its transverse energy dominates more than 0.95 of the jet. The second 
method is a more general-purpose; we define j T as a jet which contains 1 or 3 charged tracks in a 
small cone (R = 0.15) centered at the jet momentum direction with the transverse energy deposit 
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FIG. 3: Invariant-mass distributions of £ir T (left) and £j T (right) in the pp — > H ++ H~~ —> £ + £ + £~t~ process 
followed by hadronic decays of t~, after the requirement of the same-signed dilepton with M u ~ m H ±±. 
Dashed (Solid) histograms are for iTjf (H^). Smooth lines in the left panel are theoretical expectations 
by using Eqs. (j4|). 



to this small cone more than 0.95 of the jet. The second method could tag the r decay into 2tt and 
3tt originated from r —> pv and a\v decays, in addition to the single it. Thus the tagging efficiency 
is better than the first method, but the spin analysis power is weakened. 

The extensive signal-to-background studies including r's can be found, for example, in Ref. 13|, 



[jjj]. Following their results, a clear signal extraction is i 
signed dilepton and possibly a peak in their invariant-mass. Thus, we present the simulation results 
only for the signal events, but not for the background events. Expected background processes are 
diboson production, tt plus one boson production, etc. 

As the first case, we deal with pp — > H ++ H — > £ + £ + £~t~. This decay pattern can be 
easily identified by requiring the same-signed dilepton with a sharp peak in their invariant-mass 
distribution. The mass of H can be clearly obtained from the peak. Then, the invariant- 
mass distribution of the remaining £~ and decay products of r~ can be used for the polarization 
discriminant. We take the hadronic decays of r, because of the better spin analysis power than 
the leptonic ones as shown in Fig. [TJ 

In Fig. [31 the invariant-mass distributions of £tt t (left panel) and those of £j T (right panel) from 
H —> i~T~ are shown for the events with 3£ and one r-tagged jet. Results for the decay of 
(H^) are plotted in the dotted (dashed) curves. The y-axis of the plots hereafter stands for the 
number of events in our simulation. We generate the same number of events for and H^ 1 

for each decay pattern, and perform the r-tagging methods for the same event sets. Therefore, 
the difference of the total number of events between H^f 1 and H R indicates the difference of 
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FIG. 4: Invariant-mass distributions of ir T ir T (left) and j T j T (right) in the pp —> H ++ H~~ —> £ + £ + t~t~ 
process followed by hadronic decays of t - 's, after the requirement of the same-signed dilepton with M u ~ 
m H ±±. Dashed (Solid) histograms are for (H^ ). Smooth lines in the left panel are theoretical 

expectations by using Eqs. (j6|). 



efficiencies due to the kinematical cuts; the difference of the total number of events between the 
curves for ir T and j T indicates the difference of the r-tagging efficiencies between the two methods. 
However, the comparisons of the total number of events between the plots in different decay patterns 
are meaningless. The correct comparison requires the branching ratios of ff ±=t which we do not 
specify in this study. The distributions in the left panel roughly reproduce the linear dependence 
on z given in Eqs. ([4]) which is superimposed with an appropriate normalization. Thus, the pionic 
r-tagging method seems to be working well to catch the pionic r decay products. The effect of 
the kinematical cuts can appear in the small z region due to the p T cut. The discrimination of 
the two distributions may be possible with just a small number of events. When we employ the 
general r-tagging method, the expected number of events becomes 5 times larger than that for the 
case with pionic r-tagging method. Since the distributions for and H R still differ from each 
other, although the difference is weakened, this general r-tagging method works as well for our 
purpose. We note that the distributions in the right panel could be understood as the sum of the 
various hadronic r decay contributions with proper polarization dependence [ijj]. The distributions 
are scaled by the mass of H ±=t , thus it is expected that these distributions does not depend on the 
value of the mass so much. We have confirmed that quite similar distributions are obtained for the 
case with a heavier mass of H^. 

The second case is pp — > H ++ H — > £ + £ + t~t~ . We note that this process is reliable even if 
the branching ratio is too small for the LFV decay H — > £~t~ which is used in the first case. 
This collider signature has also a sharp peak in the invariant-mass distribution of the same-signed 
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dilepton, thus the signal event extraction from the background contributions would be easy. We 
analyze the case where both of r leptons decay into hadrons because of the best spin analysis power 
as shown in Fig. [2j 

In Fig. HI the invariant-mass distributions of tt t it t (left panel) and those of j T j T (right panel) 
from H — > t~t~ are shown. Results for the decay of (H^) are plotted in the dotted 

(dashed) curves. As expected, the distributions in the left panel are almost the same as Eqs. © 
even after the selection cuts. For the j T j T invariant-mass distributions, the curves for and for 
are still well separated. 

The third case is for a situation in which H — > £~£~ does not exist. Then the signature 
would be pp — > H ++ H — > £ + t + £~t~ , which seems better than pp H ++ H — -> ^ + r + r~r~. 
The momentum reconstruction of the two r's is still possible by using a collinear approximation 
for the decay of the two r's Then the mass of can be measured by the invariant-mass of 
the same-signed £r pairs. If both of the r's decay hadronically, we may suffer from the background 
contribution. In order to suppress the background contribution, we require that one r decays 
leptonically which gives the same-signed dilepton l£ T . The hadronic decay for the other r will be 
preferred for the discrimination of the polarization, although the £ + £~^£~£~ signature may be also 
exploited. 

In Fig. [5J the invariant-mass distributions of £ir T (left panel), those of £j T (middle panel) are 
shown for the event in which the momentum reconstruction is resolved by using the collinear 
approximation method. 3 The full momentum reconstruction is not necessary to obtain the plots in 
Fig. \5\ but is very effective to extract the signal events of this decay pattern. The obtained invariant- 
mass distributions of £ir T and £j T are almost the same as those in the pp — > H ++ H — > £ + £ + £~t~ 
process in Fig. O In addition, in the right panel, we plot the invariant-mass distributions of the 
same-signed tl T pair from the other side of the decays. The tl T invariant-mass distributions 
are in good agreement with Eqs. ©. Due to the good momentum resolution of leptons, the dilepton 
invariant-mass distributions could be also useful for the r polarization discrimination. 

Finally, we comment on the case where H decays only into £~ l~ or r~r~. If H decays only 
into t~£~ , it is impossible to observe any polarization phenomena from kinematical measurements. 
Then we should rely on predictions in each model about the production cross section and the 
decay branching ratios in order to distinguish models. If H decays only into r~r~ 

3 One may wonder how to select £ ~ for the collinear approximation method under the existence of £~ . One can find 
the correct one which gives the smaller difference of the reconstructed invariant-masses of the two ^ ± r ± pairs. 
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FIG. 5: Invariant-mass distributions of £ir T (left), £j T (middle) and ££ T (right) in the pp —> H ++ H~~ —> 
£ + t + £~t~ process followed by one leptonic and one hadronic decays of r's after the requirement of the 
proper momentum reconstruction by using the collincar approximation method. Dashed (Solid) histograms 
are for Hjr (H^). Smooth lines in the left and right panels are theoretical expectations by using Eqs. (01 
and Eqs. ([5]), respectively. 



gives the 4r signature. Since, there are too many sources of the missing momentum, the momentum 
reconstruction for the r's is not possible, we cannot use the invariant-mass peak at m H ±± for the 
background reduction. For the similar signatures in the context of the two Higgs doublet model, 
however, a sufficient background reduction is expected in the various channels for the r decays [3]. 
Thus, the invariant-mass distributions, for example, in the £^£^j T j T channel may be used for the 
polarization discriminant. Notice that the information on m H±± could be obtained by the endpoint 



of the invariant-mass distributions 



27] 



IV. CONCLUSIONS 

The doubly charged scalar boson appears in several new physics models, especially in the 
models to generate Majorana neutrino masses. The has characteristic Yukawa interactions 

with two left-handed charged leptons (e.g., A ±=t: from an SU(2)^ triplet field in the HTM) or two 
right-handed charged leptons (e.g., from an SU(2)/, singlet field in the ZBM) depending on 
the model. We have studied the kinematical consequences of the Yukawa interactions in order 
to discriminate these models through the determination of the chiral structure of the Yukawa 
interaction. 

At collider experiments, it is known that the polarization of the r lepton is analyzed by the 
energy fraction distributions of its decay products (vr 1 * 1 , £^, etc.). We have seen that the invariant- 
mass distributions are good analyzers of the r polarization especially for the hadronic r decays. 
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We have performed a simple Monte-Carlo simulation for decays of r's made from the decays 
of pair-produced ff ±=t . If H ++ H~~ — > £ + £ + £~t~ mode exists {£ = e,/j), the invariant-mass 
distribution of £~tt~ gives the best analysis power on r poralization. The background contribution 
can be highly reduced by requiring £ + £ + whose invariant-mass is ~ m H ±±. We have shown that 
H ++ H — > £ + £ + t~t~ mode followed by r — > txv for both r is also useful to discriminate the r 
polarization by using the distribution of the invariant-mass of 7r r 7r T . Even if there is no H ++ — > 
£+£+ , we have found that r polarization can be determined by the distribution of the invariant-mass 
Mfa r in H ++ H — > £ + t + £~t~ mode with a leptonic and a pionic r decays. The reduction of the 
background events can be achieved by requiring same-signed £i T with the collinear approximation 
method. Therefore, in these various cases for the leptonic decay of -ff ±± , we can determine the 
chiral structure of the Yukawa interaction of i^ ±± , and it will help to discriminate new physics 
models beyond the SM. 
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